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ADAPTIVE VERTEX-CENTERED FINITE VOLUME METHODS 
WITH CONVERGENCE RATES 

CHRISTOPH ERATH AND DIRK PRAETORIUS 


Abstract. We consider the vertex-centered finite volume method with first-order conform¬ 
ing ansatz functions. The adaptive mesh-refinement is driven by the local contributions of 
the weighted-residual error estimator. We prove that the adaptive algorithm leads to linear 
convergence with generically optimal algebraic rates for the error estimator and the sum 
of energy error plus data oscillations. While similar results have been derived for finite 
element methods and boundary element methods, the present work appears to be the first 
for adaptive finite volume methods, where the lack of the classical Galerkin orthogonality 
leads to new challenges. 


1. Introduction 

1.1. Finite volume method. A classical finite volume method (FVM) describes nu¬ 
merically a conservation law of an underlying model problem, which might be described by 
a partial differential equation (PDE). In particular, it naturally preserves local conservation 
of the numerical fluxes. Therefore, FVMs are well-established in the engineering community 
(fluid mechanics). Even though the FVM has a wide range of applications the numerical 
analysis is less developed than for the more prominent finite element method (FEM). There 
exist different versions of the FVM like the cell-centered FVM, which basically yields to 
a piecewise constant approximation of the unknown solution on a primal mesh. For more 
details we refer to |EGH00j . The so-called vertex-centered FVM (finite volume element 
method, box method) belongs to the other big family of FVMs, where one usually intro¬ 
duces an additional dual mesh around the nodes for the approximation. In this work, we 
focus on the lowest-order vertex-centered finite volume method (from now on only FVM) for 
some elliptic model problem in d = 2, 3. The first relevant mathematical analysis of this 
method started with the works |BR87[ IHac89( ICaiQlj . 

1.2. A posteriori error estimation and adaptive mesh-refinement. Accurate 
a posteriori error estimation and related adaptive mesh-refinement is one fundamental col¬ 
umn of modern scientific computing. On the one hand, the a posteriori error estimator allows 
to monitor whether a numerical approximation is sufficiently accurate, even though the ex¬ 
act solution is unknown. On the other hand, it allows to adapt the discretization to resolve 
possible singularities most effectively. Over the last few years, the mathematical understand¬ 
ing of adaptive mesh-refinement has matured. It has been proved that adaptive procedures 
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for the finite element method (FEM) as well as for the boundary element method (BEM) 
lead to optimal convergence behavior of the numerical scheme; see, e.g., |D5r961 IMNSnni 
IBDPndl [StiOTl KlKNSnSl IFFP14] for FEM, |FKMP131 lFFK+141 lFFK+151 Kfanl3j for BEM, 
and [CFPP14] for some general framework. 

In this work, we analyze an adaptive mesh-rehning algorithm of the type 


SOLVE 


ESTIMATE 


MARK 


REFINE 


( 1 ) 


in the frame of the FVM (Algorithm 0]). Given a conforming triangulation 7^, the module 
SOLVE uses FVM to compute a discrete approximation ui to the solution u of the PDE. For 
the ease of presentation, we assume that the linear system is solved exactly, although, in the 
spirit of |CFPP14l Section 7], stopping criteria for iterative solvers can be included into our 
analysis. The module ESTIMATE employs a weighted-residual error estimator rjg, from |CLT051 
IXZOGj which is also well-studied in the context of adaptive hnite element methods |Ste07( 
ICKNS08|IFFP14] . The module MARK uses the Dorfler marking criterion introduced in [Dor96j . 
to mark elements for rehnement, where the local error appears to be large. Unlike common 
algorithms for FEM and BEM, we follow |MNS00] and also mark elements with respect to 
the data oscillations to overcome the lack of the Galerkin orthogonality. Finally, the module 
REFINE employs newest vertex bisection (NVB) to rehne the marked elements and to generate 
a new conforming triangulation T^+i which better resolves the present singularities. 


1.3. Contributions of the present work. Iteration of the adaptive loop ([T]) provides 
a sequence of successively rehned triangulations 7^ together with the corresponding FVM 
solutions Ui and the a posteriori error estimators ri£. Theorem [3 below proves that this 
adaptive iteration leads to linear convergence in the sense of 

Vi+n < Cq'^ r]i for all £,n eNq (2) 

with some independent constants C > 0 and 0 < g < 1. Under an additional assumption on 
the marking which can be monitored a posteriori, we prove optimal convergence behavior 

Vi<C(#T,-#Tor ( 3 ) 

for each “possible” algebraic rate s > 0 (in the sense of certain nonlinear approximation 
classes which are dehned in Section [2U] below), where ^Ti denotes the number of elements 
in Ti- These results can be equivalently stated with respect to the sum of energy error plus 
data oscillations, which is usually done in the FEM literature |Ste07( IGKNSOSl IFFP14j , since 

C~^rie < min (|||m-U£||| -Losc^(u£)) < |||m-M£||| + osc^ini) < Crjp, (4) 

•Vi 

see Theorem [2] below. We note that (|3|) in particular provides a generalized Gea lemma 
which states that the FVM solution ui is quasi-optimal with respect to the so-called total 
error, i.e., the sum of energy error plus data oscillations. Since (jl]) is also known for the 
FEM (see, e.g., |FFP141 Lemma 5.1]), this reveals that FEM and FVM lead to equivalent 
errors in the sense of 

~ '*^^111 + osc£(m^)) < |||m — -L osc£(m^™) < C (|||m — Ui\\\ + oscf(M£)), (5) 

where is the FEM solution with respect to the FVM space. This complements recent 

results which compare the total errors of different FEM discretizations |GPS12l IGKPS15] . 
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Unlike the results for FEM and BEM, the novel Cea-type estimate (jl]) as well as our 
result ©-([S]) on adaptive FVM requires the additional assumption that the initial triangu¬ 
lation To is sufficiently fine. We note, however, that such an assumption is also required to 
prove well-posedness of the FVM in general and thus appears naturally. 

Prior to this work, a posteriori error estimates for the FVM for elliptic model problems are 
derived in |CLT051 fXZOGl IZoulOj : see also |Eral31 Remark 6.1] and |Eral31 Conclusions] for 
estimates which are robust with respect to the lower-order convection and reaction terms. 
To the best of the authors’ knowledge, convergence of an adaptive 2D FVM has only been 
analyzed in the yet unpublished preprint |XZ06j . The latter is concerned with convergence 
only and the analysis follows |MNS00] and relies on a discrete efficiency estimate and hence 
on the so-called interior node property of the mesh-refinement. Contrary to |XZ06] . our 
analysis extends the ideas of |CKNS08] and provides a contraction property for the weighted 
sum of energy error, weighted-residual error estimator, and data oscillations. Therefore, 
our analysis covers in particular standard NVB, where marked elements are refined by one 
bisection. 

We hnally note that residual error estimators have also been developed for the cell-centered 
hnite volume method |Nic051 lEPOSt IVohOS] . These a posteriori estimators rely on an inter- 
polatory post-processing of the original piecewise constant cell-centered hnite volume ap¬ 
proximation. Thus, a thorough adaptive convergence analysis requires additional ideas to 
extend and adapt the analysis presented below. 

1.4. General notation. We use < to abbreviate < up to some (generic) multiplicative 
constant which is clear from the context. Moreover, ~ abbreviates that both estimates < 
and > hold. Throughout, the mesh-dependence of (discrete) quantities is explicitly stated 
by use of appropriate indices, e.g., Ux is the FVM solution for the triangulation Tx and 77 ^ 
is the error estimator with respect to the triangulation Ti. 

2. Model problem & main results 

2.1. Model problem. Let D C M'^, d = 2,3, be a bounded and connected Lipschitz 
domain with boundary F := dVL . As model problem, we consider the following stationary 
diffusion problem: Given / G LF‘{VL), End u G H^{VL) such that 

—divAVw = / in D, (6a) 

M = 0 on F. (6b) 

We suppose that the diffusion matrix A = A(a;) G is bounded, symmetric, and uni¬ 
formly positive dehnite, i.e., there exist constants Amin, Amax > 0 such that 

Amin |vp < v^A(a;)v < Amax |vp for all V G M'’* and almost all x G D. (7) 

For convergence of our FVM and well-posedness of the residual error estimator, we addition¬ 
ally require that A(x) is piecewise Lipschitz continuous, i.e., 

A G for all T G To, (8) 

where To is some given initial triangulation of D; see Section [231 below. 

The weak formulation of the model problem ([6]) reads: Find u G Hq{Q) such that 

A{u,v) := (AVm, Vt)q = if ,v)n for all v G 
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where (0,'0)n = dx denotes the L^-scalar product. According to our assump¬ 

tions o on A, the bilinear form A(-,-) is continuous and elliptic on Hq{Q). Therefore, 
existence and uniqueness of the solution u G Hq{Q) of ([9]) follow from the Lax-Milgram the¬ 
orem. Moreover, |||f||p := A{v,v) defines the so-called energy norm which is an equivalent 
norm on HliVl). We shall use the notation 111^111^ := f^AVv ■ Vv for the energy norm on 
subdomains a; C hi, i.e., |||n||| = |||w|||o. According to ([7]), it holds |||n|||w — ||Vn||L 2 (^^). 


2.2. Triangulation. Throughout, 77 denotes a conforming triangulation of hi into non¬ 
degenerated closed simplices T G 77 (i-e., triangles for d = 2, tetrahedra for d = 3), Afx 
is the corresponding set of nodes, and J^x is the corresponding set of facets (i.e., edges for 
d = 2 and triangular faces for d = 3). We suppose that 77 is a-shape regular, i.e.. 


diam(T)'^ 

max-—- < a < oo. 

T€Tx |T| 


( 10 ) 


Here, diam(T) ;= max {|a;—1/| : x,y E Tj denotes the Euclidean diameter and |T| is the area 
of T. Additionally, we assume that the triangulation 77 is aligned with the discontinuities 
of the coefficient matrix A, i.e., ([H]) holds with % replaced by 77- We note that this follows 
from dH]) and the mesh-refinement used; see Section 12.61 Associated with 77 is the local 
mesh-size function hx G L°^(Q) which is dehned by hxlr ■= hx ■= Note that a-shape 

regularity (ITU]) yields hx — diam(T). 

For the nodes Afx , we introduce the partition A/7 = AT^ U Af^ into all boundary nodes 
A^x := Afx n r and all interior nodes Af^} := Afx\Afx- 

For the facets J-'x, we introduce the partition A'x = J^x ^ boundary facets 

A'x •= ^ : A C r} and all interior facets := J^x\J^x- Finally, for an element 

T G 77, we denote by AV := {F G Fx '■ F C dT } C Fx the set of all facets of T. 


2.3. Dual mesh. In contrast to standard FEM, our FVM discretization additionally 
needs the so-called dual mesh Ff which is built from the conforming triangulation 77- In 
2D, connecting the center of gravity of an element T G 77 with the (edge) midpoint of 
F G Ft, we obtain Tf whose boxes (elements) V E Tf are non-degenerate closed polygons; 


see Figure [11(a) In 3D, we connect the center of gravity of an element T G 77 with the 
centers of gravity of the four faces F E Ft- Furthermore, each center of gravity of a face 
F G Ft is connected by straight lines to the midpoints of the edges of the face F. Figure [ ^aj] 
shows the contribution of some element T eTx with node Oj to the box Vi E Tf. 

Note that there is a unique correspondence between the nodes a, G Afx of the primal mesh 
77 and the boxes Vi E Tf of the dual mesh, namely 17 n ATx = {oj}. Furthermore, we define 
;= {F n 1/ : F G J^x} for all F G 77; see Figure [fb)] for 2D. 

For 3D, Figure l ^b)] shows three faces Ci, C 2 , and (3 of Fvi,x, Vi E 77- Note that 


U T = Sl= U K and U f’= U U (11) 

T&Tx V&T* F&Tx VgT* F€Fv,x 


2.4. Vertex-centered finite volume method (FVM). Given the conforming trian¬ 
gulation 77 and the corresponding dual mesh 77, we define the space of all 77-piecewise 
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Figure 1. Local construction of the dual mesh 7^ from the primal mesh 
7^ in 2D: The dashed lines are the boundaries of the induced control volumes 
Vi G 7^, which are associated with the nodes Oj G A/”x of Tx (left). For V G 
the set J^v,x consists of the bold lines which are parts of edges in (right). 

affine and globally continuous functions 

<S^{Tx) ■= {n G C(r2) : v\t affine for all T G 7)<} C 77^(17) 

as well as the space of all T^-piecewise constant functions 

:= {n G L^iyt) : v\v constant for all F G 7^}. 

For the FVM discretization, we consider the subspaces which respect the homogeneous 
Dirichlet conditions of (j6]), i.e., 

Sl{Tx) := {v G S\Tx) : i;|r = O} C Hl{n) and P°(r;) := {v G V\T:) : v\r = O}. 

The formal idea of the FVM reads as follows: If we integrate the strong form ([6]) over each 
dual element V and apply the divergence theorem, we get a balance equation for the 
model problem. The FVM approximates u G 77o(D) by some conforming approximation 
Ux E SliTx) of the balance equation. With the aid of test functions in Vq{Tx)-i we formalize 
this with the bilinear form 

Ax{vx,v*^) := - V] v*^\vi [ AVux ■ nds for all Ux e Sq{Tx) and n* G Vq{T*). qo! 
The right-hand side reads 

<lu / fdx = {f,v*Jn for all < G Po°('^x)- 

Throughout, if n appears in a boundary integral, it denotes the unit normal vector to the 
boundary pointing outward the respective domain. Now, the FVM discretization reads: Find 
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Figure 2. Local construction of the dual mesh 7^ from the primal mesh 7^ 
in 3D: For a node a, G Afx of T, the center of gravity ct of T is connected 
with the centers of gravity cp. of the three adjacent faces F G Ft- Moreover, 
these centers are connected to the midpoints of the three edges which meet 
in Oj. Together with the edges from these midpoints to a*, we get the cuboid 
Vi DT 7 ^ 0 (left). The three dark-gray faces Qi, Q 2 , and Q 3 are part of the 
boundary dVi of the box Vi (left). The light-gray faces Ci, C 2 , and (3 belong to 
the set Fvi,x and are part of faces in Fx (right). 


Ux ^ Sq(Tx) such that 

^x{ux,v*J = {f ,v*^)n for all n; G T’o(Tx). (13) 

It is well-known that there exists a constant H > 0 such that flT^ admits a unique solution 
Ux G Sq{Tx) provided that Tx is sufficiently hne, i.e., ||hx||Lo°(o) < H] see Lemma HH below. 
The convergence of the FVM is usually proved under certain regularity assumptions, e.g., 
u G H^{n) n H^+%n) for some e > 0 ; see, e.g., |ELL02( Theorem 3.3.]. As a side result of 
our analysis. Theorem [3] below proves convergence of the total error (i.e., energy error plus 
data oscillations) without any regularity assumptions. 

2.5. Weighted-residual error estimator. With divx denoting the 7^-piecewise di¬ 
vergence operator, we dehne the volume residual by 

-Rx(nx)|r = (/+ divxAVnx)|r for all T G Tx and all Ux e 5o(Tx). (14) 

Throughout, we abbreviate divx AVux := divx (AVux) to ease the readability. Let [■] denote 
the normal jump across an interior facet T = T n T G F^, i.e., [^] |f = 9\t ■ n-r + qIt' ■ t^t', 
where, e.g., qIt denotes the trace of g from T onto F and uj’ is the outer normal of T on 
F. Then, we dehne the facet residual or normal jump by 

F{vx)\f = [AVuxIIf for all F e F^ and all Ux e 5o^(rx). 
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For all Vx G iSq( 71<), we define the weighted-residual error estimator as for the FEM 

?7x(wx)^ = ?7x(Tx,nx)^ with ^ forallWx^Tx, (16) 

TeUx 

where 

Vx{T,Vx) = h-p 11/-|-divxAVux 11^2(2’) + hr II[AVux]||L2(g2’yp); (17) 

cf., e.g., |AOOOl IVerl3] . For Vx = Ux being the discrete FVM solution, we abbreviate the 
notation and omit this argument, e.g., rjx ■= rjxiux) and rjxiT) := ri(T,Ux)- 
Let fix denote the elementwise or facetwise integral mean operator, i.e.. 



Recall that fix is the elementwise L^-orthogonal projection onto the constants, i.e., 

||n — nxn||i, 2 (T-) = min ||n — c||i2(.r) ^ II'^||l 2 (t-) for all r G 7j< U J^x and all v G (19) 

With fix, we dehne the data oscillations 

oscx(nx)^ = osCx(Tx,nx)^ with osCx(WxWx)^ = ^ osCx(T,Ux)^ for all Ux ^Tx, 

T&Ax 

where 

osCx(T,nx)^ = 11(1 - nx)(/ -h divxAVnx)||i2(r) -h hr ||(1 - nx)[AVnxllli2(ar\p). (20) 

Again, we abbreviate the notation for Vx = Ux being the FVM solution, e.g., oscx := 
oscx(nx) and osCx(T) := osCx(T, Ux)- Moreover, we stress the elementwise estimate 

osCx(T,Ux) < Vx{T,Vx) for all T G 7)< and all Vx G >So(Tx) (21) 

which immediately follows from (IT^ . The following result is proved in |XZ06l Theorem 2.4] 
and |XZ061 Theorem 2.6]; see also |CLT05l Theorem 3.1] and |CLT05l Theorem 3.3]. 


Proposition 1 (reliability and efficiency). The estimator rix satisfies reliability 

lll„, ll|2 , 


as well as efficiency 


^eflf Vx 


'reivl 

(22) 

p OSC^. 

(23) 


The constants Crei, Cefr > 0 depend only on a-shape regularity of Tx and on the assump¬ 
tions (I7])-(IH]) on A. □ 


The hrst contribution of the present work is the following Cea-type quasi-optimality of 
FVM with respect to the total error (i.e., sum of energy error plus data oscillations). In 
particular, this implies that the total errors of FVM and FEM are equivalent, see (]5]). The 
proof of the theorem is given in Section 13.61 


Theorem 2. There exists H > 0 such that the following statement is valid provided that Tx 
is sufficiently fine, i.e., ||hx||L°°(o) < H: There is a constant Ctot > 0 such that 

Cfi,lVx< min (|||n - Ux III + osCx(ux)) < |||m - Mx III + osCx < Ctot ?7x- (24) 

vxes^(Tx) 
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Figure 3. 2D newest vertex bisection: For each triangle T G T, there is one 
hxed reference edge, indicated by the double line (left, top). Rehnenient of 
T is done by bisecting the reference edge, where its midpoint becomes a new 
node. The reference edges of the son triangles are opposite to this newest 
vertex (left, bottom). To avoid hanging nodes, one proceeds as follows: We 
assume that certain edges of T, but at least the reference edge, are marked for 
rehnement (top). Using iterated newest vertex bisection, the element is then 
split into 2, 3, or 4 son triangles (bottom). 


Moreover, if 6 iSq (7)<) denotes the FEM solution of 

^(Mx™,Wx) = {f,Vx)n for allvx G <So{Tx), (25) 

it holds 

Qot (Ilk -“x III +OSCx) < ||k-Mx™||| +OSCx(Mx™) < Ctot (||k-Mx||| +OSCx). (26) 

The constant Utot > 0 depends only on D, H, the a-shape regularity of Tx, and on the 
assumptions dZD-® on A. 

For the sake of completeness and as an application of Theorem [2], we note the following 
a priori estimate for the total error. Note that fl27)) does not require any additional regularity 
assumption on u. The proof is given in Section 13.71 


Theorem 3. There exists H > 0 such that the following statement is valid provided that Tx 
is sufficiently fine, i.e., ||/ix|| l°°(o) < H: There is a constant C > 0 such that 

U'^dllM-Mxlll+osCx) < ||hx(l-nx)/||L2(t^) + min (||k - Wx III + IkxVux ||l2(o))• ( 27 ) 

Vx&S^{Tx) 


In particular, this proves convergence 

Ilk — WxIII + oscx 0 as Ikx|| l°°(o) 0. (28) 

Provided that u G Hq{Q) H there even holds 

|||m-Mx||| +oscx = C’(|kx||L-(o))- (29) 

The constant C > 0 depends only on H, H, the a-shape regularity ofTx, and on the as¬ 
sumptions o -([8]) on A, and fl28|) - fl2^ require uniform a-shape regularity of the considered 
family Tx ■ 


2.6. Adaptive algorithm &; main result. As for adaptive hnite element meth¬ 
ods |Dor96l IMNSOOl ISte071 ICKNS081 IFFPldj . we consider the following adaptive algorithm 
which specihes the adaptive loop ([1]). Unlike the common algorithms in the context of adap¬ 
tive FEM and BEM |CFPP14] . our algorithm does not only employ Dorfler marking with 
respect to the error indicators rjiiT), but also for the local contributions osci(T) of the data 
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oscillations. This additional marking step is necessary to control the lack of Galerkin orthog¬ 
onality (138|) and thns allows to prove (linear) convergence fl3^ of the adaptive algorithm. 

For the mesh-rehnement in step (v) of Algorithm 01 we employ newest vertex bisection 
(NVB); see, e.g., |KPP131 ISteOSj for general dimension d>2 and Figure 0] for an illustration 
for d = 2. For a conforming triangulation T and a set of marked elements AT C 7~, let 
T' := refine(T, AT) be the coarsest conforming triangulation generated by NVB such that 
all marked elements T G AT have been rehned, i.e., AT C 'T\T'. 

Algorithm 4. Input: Let 0 < 6*' < 6* < 1 and Gmark > 1 given adaptivity parameters. 
Let To be a conforming triangulation of VL which resolves possible discontinuities of A in the 
sense of (|8]). 

Then: For i = 0,1,2,... iterate the following steps (i)-(v).' 

(i) Solve (0^ to compute the discrete solution U£ G Sq{T) corresponding to %■ 

(ii) Compute the refinement indicators rjfiT) from (ITTI) and the data oscillations from fl20|) 
for all T E T- 

(hi) Construct a subset AT^ ^ T of up to the multiplicative factor Gmar k minimal cardi¬ 
nality which satisfies the Ddrfler marking criterion 

evI<Vi{M^£Y- (30) 

(iv) Construct a subset Aii F % of up to the multiplicative factor G^ark minimal cardi¬ 
nality which satisfies Ai] C AT^ as well as the Ddrfler marking criterion 

9' oscj < oscfiMifi. (31) 

(v) Define T+i := ref ine(7£, AT^). 

Output: Adaptively refined triangulations Ti, corresponding discrete solutions U£, estimators 
Tji, and data oscillations osc^ for £ > 0. 


Remark 5. (i) For Gmar k = 1, the construction of the set ATI in step (iii) of Algorithm\^re- 
guires to sort the error indicators and thus results in logarithmic-linear complexity. Instead, 
for G mar k = 2 , ou approximuto sorting based on binning allows to construct Aic in linear 
complexity |Ste07j . The same applies for At £ in step (iv) of Algorithm\^ 

(ii) There exists a constant H > 0 such that (na has a unigue solution provided that 
\\h 


|l°°(o) < H; see Lemma\^below. Since NVB guarantees ||h£|| 2 ,oo(Q) < ||ho||L<x>(r 2 ), it 
is sufficient to suppose that the initial triangulation To is sufficiently fine. 

(iii) In step (v) of Algorithm one may use any variant of NVB which applies at most n 
bisections per marked element, where n > 1 is a fix constant. □ 


Next, we dehne certain nonlinear approximation classes, which are needed to prove optimal 
convergence behavior ([3]). To this end, we write T' G refine(T), if there exists some 
n G No, triangulations T-, and marked elements AT) C T- such that T = Tf, T' = Tf, and 
Tj = ref ine(7~'_i, AT'_^) for all j = l,...,n. Given To from Algorithm 01 we note that 
NVB ensures that all triangulations Tx G refine (To) are uniformly a-shape regular ffTOj) . 
where a depends only on Tq. 

For A > 0, we abbreviate Tat ;= [Tx ^ refine(To) ; fiTx — fiT < A}, where ft^Tx 
denotes the number of elements in Tx- For all s > 0, we dehne the approximability measure 

||m|1a, := sup min {N + Viffix, 

N>0 TxGTjv 
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where rjx denotes the weighted-residual error estimator ffT^ associated with the optimal 
triangulation Tx- Note that ||m||as < oo means that an algebraic decay rjx = 0{N~^) is 
theoretically possible if for each iV > 0 the optimal triangulations 7^ G T^r are chosen. 

As a corollary of Theorem |21 we obtain that the corresponding approximation class (of 
all u which satisfy ||m|1as < oo) can equivalently be characterized by the so-called total error 
(i.e., energy error plus data oscillations) and hence coincides with the approximation classes 
from the FEM literature; see, e.g., |Ste07l ICKNS081 IFFPldj . 

Corollary 6. There exists H > 0 such that the following equivalence is valid if the initial 
triangulation % is sufficiently fine, i.e., ||ho||i,oo(f 7 ) < H: For all s > 0, it holds 

||m||as < C )0 sup min inf (A^-|-1 )®(|||m — nx||| + oscx(nx)) < C) 0 . 

V>oT<GTjv Vx€S^{Tx) 

Proof. Note that all triangulations Tx G refine (To) are uniformly cr-shape regular and 
satisfy ||hx||Loo(Q) < ||ho||L°o(o). Therefore, the claim follows from fl2T)) . □ 

Besides Theorem [21 the following theorem is the main result (|2])-([3]) of our work. Un¬ 
like |Ste07[ ICKNS08] . we follow |CFPP14l IFFPldj and formulate the result with respect to 
the error estimator as this is the natural goal quantity of Algorithm jH In view of 0241) . the 
theorem can equivalently be formulated with respect to the total error. Its proof is given in 
Section 13.91 below. 


Theorem 7. There is a constant H > 0 such that the following statements (i)-(ii) are valid 
provided that the initial triangulation To is sufficiently fine, i.e., ||ho||Loo(Q) < H: 

(i) For all 0 < 9’ < 9 < 1, there exist constants Cun > 0 and 0 < gun < 1 such that the 
adaptive Algorithm^ guarantees linear convergence of the estimator in the sense of 

Ve+n < C'linTun Ve for all i,n e No- (32) 

(ii) There exists a bound 0 < 9opt < 1 such that for all 0 < 9 < 9opt, the following holds: 
Provided that there is a constant Cmns > 1 such that fiM.i < Cmns#A4^ for all 
i G No, there is a constant Copt > 0 such that for all s > 0, it holds 

^l+s 

||m||a, < OO \\u\\ks {f^Ti - (33) 

(1 - Tlin )* 

i.e., the adaptive algorithm leads asymptotically to each possible algebraic decay s > 0 
of the error estimator. 

The constant 6'opt depends only on hi, Ft, and uniform a-shape regularity of the triangulations 
Tx G refine (To), the constants Cun and gun depend additionally on 9 and 9', while the 
constant Copt depends also on the use of NVB and on Cm ar k and Cmns- 


Remark 8. (i) The additional assumption in Theorem 0(i>) assumes that marking ([3l]) of 
the data oscillations is negligible with respect to the overall number of marked elements. We 
note that 9' > 0 can be chosen arbitrarily small so that, in practice, fl30|) already implies dll. 
(ii) Instead of the additional marking step (iv) in Algorithm\^ one can also define M-a := 
and monitor a posteriori if 


. t oscfiMifi , 
sup ml - 5 - =: 9 > 0. 

£oGNo ^^^0 OScj 


( 34 ) 
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In this case, linear eonvergenee fl3^ with optimal rates (l3^ follows. However, for O' = 
0, even convergence remains mathematically open, so that we favor the present form of 
Algorithm\^ which guarantees (l32|) . while fl33l) requires an additional assumption. □ 

3. Proofs 

3.1. Axioms of adaptivity. In [CFPP141 Theorem 4.1], it is proved in a general 
framework that the following set of four axioms is sufficient (and partially even necessary) 
to guarantee linear convergence with optimal algebraic rates in the sense of Theorem [3 In 
particular, the model problem, the discretisation, and the estimator enter only through the 
proof of these axioms. Implicitly, we assume that given Tk G refine (To), the corresponding 
FVM solution Uk G S^lTk) is well-defined. With this convention, the axioms read: 

(Al) stability on non-refined elements: There exists a constant C > 0 such that for 
all To G refine (To) and all Tx G refine(T^), it holds 

\vx{Tx n%) -rjoiTx nT;)| < C|||mx -MoIII- 

(A2) reduction on refined elements: There exist constants 0 < g < 1 and C > 0 such 
that for all T G refine (To) and all Tx G refine(T^), it holds 

Vx{Tx\Tf < qy]o{%\Txf + c |||mx - Mo|||^. 

(A3) general quasi-orthogonality: There exists C > 0 such that for all £ G No, h holds 

OO 

\\\uk+i - Uk\\\^ < Crjj. 
k=e 

(A4) discrete reliability: There exists a constant C > 0 such that for all T G ref ine(To) 
and all Tx G refine(T^), there exists some set IZo C % with %\Tx C To and 

ifTo < Cif{To\Tx) as well as |||mx - tto||| < Crj^iTo). 

The subsequent analysis proves that Algorithm 0] for our adaptive FVM guarantees the 
validity of (A1)-(A4) if the initial triangulation To is sufficiently fine. 

3.2. Stability reduction of error estimator. The following lemma is stated without 
a proof, since the details are implicitly found in |CKNS08l Section 3.1]. Moreover, (Al)-(A2) 
do not only hold for the FVM solutions Ux and Uo (or for FEM solutions in [CKNS08] h but 
for all Vx G SliTx) and Vo G iSo(T^). The reader may simplify the proof of Lemma fTOl below, 
which provides slightly sharper estimates for the data oscillations. 

Lemma 9. The residual error estimator satisfies the following two properties (Al’)-(A2’).- 

(AP) There exists a eonstantC > 0 such that for allTo G refine (To), allTx G refine(T^), 
and all G Sq{Tx), Vo G iSq(T^), it holds 

1 /2 

\vxiTx nTo,vx) - voiTx nTo,vo)\ < C lll^x - woIIIt) • 

TG^x {'Tx riTo) 


11 








(A2’) There exist constants 0 < q < 1 and C > 0 such that for all To G refine (To), all 
Tx G refine(7^), and all Vx G Sq{Tx), Vo G Sq{To), it holds 

VxiTx\T,Vy,)‘^ < qr]o{T\Tx,Vof + C ^ \\\vx - Vo\\\l^ ^ . 

T€nx{Tx\%) 

Here, VLx{Ux) '■= {T &Tx '■ TT' G Wx TflT' 7 ^ 0} denotes the patch ofUx ^ Tx inTx- 
The constants C and q depend only on uniform a-shape regularity of all Tx G ref ine(7o) 
and the assumptions ([7|) -l|5 |) on A. In particular, this implies (Al)-(A2). □ 

3.3. Stability Sz reduction of data oscillations. Our proof of linear convergence fl32|) 
in Section 13.41 requires to control the data oscillations which arise in some quasi-Galerkin 
orthogonality 038 p . This is done by means of two additional axioms which structurally fol¬ 
low (Al’)-(A2’), but have an additional factor ||hx||Loo(o) in the perturbation term. Es¬ 
sentially, the following lemma is a sharper variant of the proofs in |XZ06[ Lemma 5.2] 
and |CKNS08l Section 3.1]: 

Lemma 10. The data oscillations satisfy the following two properties (B1’)-(B2’).- 

(Bl’) There exists a constant C > 0 such that for allT G refine (To), allTx G refine (7^), 
and all Vx G iSq(7^), G iSq(7^), it holds 

1 y/2 

|oscx(Tx n 7^,nx) - osco(Tx n 7^,no)| < G ^ ^ -'^<>111^) 

(Tx nTo) 

(B2’) There exist constants 0 < q < 1 and C > 0 such that for all T G refine (To), all 
Tx G refine(7^), and all Vx G Sq{Tx), Uo G Sq{T), it holds 

oscxiTxXT^^Vx)"^ < qosco{T\Tx,VoY + C ^ ^rlll^x - WoIIIt) "" 

T(^nx{Tx\%) 

Here, Qx{h(x) '■= {T & Tx ■ 3T' G Wx T n T' 7 ^ 0} denotes the patch of Ux ^ Tx in 
Tx- The constants C and q depend only on uniform a-shape regularity of the triangulations 
Tx G refine (To) and the assumptions (Tfp -llSP on A. 

Proof. Step 1. For all Wx G iSo(T^) and all Ux ^ T^, it holds 

h^||(l - nx)divxAVwx|li2(r) < G ^Tl|Vwx|li2(r), (35) 

TgWx TOUx 

where G > 0 depends only on maxTijgT^ || A IIFor T G Wx, it holds 

hrlKl — hlx) divxAVwx||L2(r) < ^t|| divxAV wx||l2(t) ^ ^r||A||vi/i.°°(T)||Vr(;x||L2('r), 

since Vwx is constant on T. All elements T E Tx satisfy T Tq for some Tq G To, i.e., 
||A ||wi.o°(r) < max-ToGT;, ||A|| Summing this estimate over all T G Ux, we thus 
obtain ([35]). 

Step 2. For all Wx G SliTx) and all Ux ^ Tl<, it holds 

5 ; ;,rll(i-n 0 Jx ('^x)||i2(g7’\r) — ^ E /^rll Vwx |Il 2 ( 2 -), (36) 

T&Ax Tenx(Ux) 
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where C > 0 depends only on m.ax.To&To II A|| vyr°°(To) and a-shape regularity of Tx '■ Let T &Ux 
and F G Ft H F^ be a facet of T which is not on the boundary L. Let A = (1/|T|) A dx, 
i.e., piecewise integral means of the entries in A. Note that VtCx as well as the outer normal 
vector ht of T are constant on F. The uniform continuity of A It, the Poincare inequality 
in and a scaling argument show 

- nx)(AVM;x ■ nT)||L2(T) < h^'^^IKA - A)Vm;x ■ nT||L2(T) 

< IIA — A||too(t) hj; ||Vr(;x||L2(F) 

< hr IIA||vi/i>°°(t) II Vu’x||l2(t). 

Let T' G 7^ be the unique element with F = T (1 T'. Then, the definition of the facet 
residual flTS]) on F leads to 

hrlKl — Llx) Jx('M^x)||i2(T) ^ h-T || Vwx ||l2(t) + ^r' ll^'^x ||i2(T')- 

Summing this over all interior facets F G Ft H F^ of the elements T eUx, obtain (136|) . 
Step 3. For all Vx,v'^ G iSq( 7^) and all Wx C 7^, it holds 

|0SCx(Wx,Ux) - OSCx(Wx,Ux)l < ^ II “ l^x) II L2(T) j > (37) 

TeS^x {Ux ) 

where C > 0 depends only on maxTocTs ||A||v[/i. 00 (^ 0 ) and a-shape regularity of Tx- The 
inverse triangle inequality for square-summable sequences in the Banach space £2 gives 

|oSCx(Wx,Ux) - OSCx(Wx,t’x)| <(^ ^ 11(1 - fix) diVxAV(Ux - Wx)||i2(T) 

T&Ux 

\ 1/2 

+ / , 11(1 ~ Ifx) Jx(llx ~ llx)llL2(9T\r)j 

T&Ax 

Using fl35ll - fl^ for Wx '■= Vx — v'^, we obtain ([37]). 

Step 4: Proof of (Bl’). For Vx G iSq( 7^) and Vo G iSo(77), apply fl37)) with v'^ := Vo 
and Ux '■=Tx^T- Note that osCo(Tx n7^, Uo) = osCx(7^ n77, t'o)- With || V ■ ||i;, 2 (T) ~ ||| ■ |||t 
this yields 

1 

|oSCx(Tx n7;,Ux) -OSCo(Tx n7;,Uo)| < (^ ^ ^tIII^x - IIoIIIt) • 

TeS^x (7x nTo) 

The hidden constant depends only on maxToGT^ ||A||vi/i,oo(Tq), on a-shape regularity of Tx, 
and on the assumptions (I7])-([8]) on A. This concludes the proof of (Bl’). 

Step 5: Proof of (B2’). For Vx € iSo(77<) and Uo £ *50(77), we apply (j37]) with v'^ := 
and Ux '■= Tx\%- With || V • ||l 2 (t) — ||| • |||t this shows 

1 y/2 

OSCx(Tx\7;,Ux) < OSCx(Tx\7;,Uo) + E ^ ^tIII^Ix - Wollll) • 

TGOx(Tx\ro) 

For all (5 > 0, the Young inequality (a + 6)^ < (1 + 5) + (1 + 6 ~^) b‘^ proves 

oscx(Tx\7;,Ux)^ < (l + (5)oscx(Tx\7;,Uo)^ + (l + 5"^)E^ ^ /i||||ux - WoIUt- 

TGOx(rx\ro) 
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Let T G 7l\Tx- Let Tx\t ■= {T' E Tx : T' ^ T} be the set of its successors in Tx- Let 
T' e Tx\t- Recall that bisection ensures |T'| < |T|/2. With 0 < g := < 1, it follows 

oscx {T , Uo) < q (^hrp \\(1 — no)(/ + divoAVno)|| l 2 (' 7 ’/) + hx ||(1 — 11^)[AVuo]||L 2 ((g 7 ’/pig 2 ’)yp)), 
since AVuo is smooth inside of T so that all normal jumps inside of T vanish. This leads to 

osCx(Tx\7;,no)^ = ^ osCx(T',no)^ < g ^ osCc{T,v^f = qosCo{%\Tx,Vof. 

TeToVTx T'sTxIt TsToVTx 

Choosing 5 > 0 such that 0 < g := (1 + 5) g < 1, we conclude the proof of (B2’). □ 

3.4. General quasi-orthogonality &: linear convergence. The following proposition 
proves (1521) in TheoremI7|(i) and shows, in particular, that the general quasi-orthogonality (A3) 
is satished. 

Proposition 11. There is a constant H > 0 such that the following statement is valid pro¬ 
vided that ||ho||L°o(Q) < H: For all 0 < 9' < 9 < 1 Algorithm^ guarantees linear convergence 
in the sense of Theorem I71(i). Moreover, together with reliability (1221) . estimate (1521) also 
implies the general quasi-orthogonality (A3). 

Our proof relies on the following quasi-Galerkin orthogonality property from [XZ06] . 
Lemma 12. Let Tx E refine(7^). Then, the corresponding discrete solutions satisfy 

|||m — Ux IIP < |||n — M^IlP “ (1 ~ ^) lll'^x — w^IlP + ^'^GgaiosCx (38) 

for all 0 < 6 < 1. The constant Ggai > 0 depends only on a-shape regularity ofTx ■ 

Proof. According to |XZ06l Theorem 5.1], it holds 

\A{u — Ux,Vx)\ < c llluxlll osCx for all Vx E Sq{Tx), 

where G > 0 depends only on a-shape regularity of Tx- For each (5 > 0, the symmetry of 
A(-, •), the last estimate, and the Young inequality 2ab < 5of -|- yield 

|||m — WxIlP = ll|■*^~■*^^llP~2 A{u — Ux,Ux — ui) — lllwx “"^^^llP 

< |||m — n^IlP — (1 — 5) |||mx — "^^^llP + osc^- 

This concludes the proof with Ggai = G^. □ 

Proof of Proposition 1771 Step 1. There exist constants Gest > 0 and 0 < gest < 1 which 
depend only on 0 < 6* < 1 and the constants in (Al)-(A2), such that 

vj+i < gest Ve + C'est \\\ue+i - n^llp for all £ eNq: (39) 

The combination of (Al)-(A2) yields for all e > 0 that 

vl+i ^ Q + (1 + Vii'Ti n T^+i)^ -|- (G -|- (1 -|- £ ^)G^) llluf+i — w^Hp. 

Note that rj{Ai]) < rp{T(\T+i). Therefore, the Dorfler marking (1501) yields 

-\-e)r]i{J'ir\T(,+\f‘ = {1 -E s)pl — {1-\-e — q)r]i{J'(\TiJ^iY 

ll30t 

< {I + e - 9 {1 +e - q))ril. 

For sufhciently small £ > 0, we see 0 < gest '■= 1 + e — 9 {1 s — q) < 1 and conclude (1501) . 
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Step 2. There exist constants Cest > 0 and 0 < gest < 1 which depend only on 0 < 0' < 1 
and the constants in (B1’)-(B2’), snch that 

osc^_^i < gestosc^ + Cest ||h^+11| ^oo 11 |m£+i - Ui\\\‘^ for all £ e Mo : (40) 

The proof follows verbatim to that of 0391) . bnt now involves (B1’)-(B2’) in combination 
with the Dorfler marking 0311) for the data oscillations. 

Step 3. Without loss of generality, we may assume that the constants Cest > 0 and 
0 < gest < 1 in fl39|) - fl40l) are the same. With free parameters 7 , /i > 0 which are fixed later, 
we define 


Av := IIIm — Us 


r + 7??x T/iosc^ • 


We claim that there are constants 7, p, C > 0 and 0 < gun < 1 such that 

< giinA^ — (1/4 — C ||h£+i||^oo(Q)) |||m£+i — MfllP, 


(41) 


where 7,/i,C, gun depend only on 6 , 9', uniform a-shape regularity of the triangulations 
Tx e refine (To), and the assumptions ([7])-([HD on A: To prove this, we use the quasi- 
Galerkin orthogonality fl55]) with 5 = 1/2 and the estimates fl59]) - (l4UD to see 


A^+I < |||m - + 77I+1 + (/i + 2Cgai)osc^+i - (1/2) |||Mf+i - 


< |||m - + gest7% + Test 


p + 2Cj 


gal 


/i 


/iOSC^ 


■( 1/2 7^est (h-+ 2Cgal)Cest ||h£+i||j;^oo(Q)) IIIm^+i M£||| 


For all £ > 0, reliability (122|) implies 

|||m - + gest 7h£ < (1 - £) |||m - Ul\\\^ + (gest + 7"^C'reie) 7T£- 

We choose 7 > 0 sufficiently small such that yCgst < 1/4. Additionally, we choose e > 0 
sufficiently small and /i > 0 sufficiently large such that 

0 < gi := gest + 7"^C'reC < 1 and 0 < g2 := gest ^ < 1. 

h 

Combining the latter estimates with C := (/i + 2Cgai)Cest > 0, we arrive at 

A^+i < (1 -e) \\\u-Ui\\f + qijr]‘^ + g2/iosc^ - (1/4 - C ||h£+i||^oc(f2)) \\\ue+i - Ui\\f 

< giin A^ - (1/4 - C ||h^+i|lL°°(o)) Ill'^^+I ~ ■*^^111^5 

where 0 < gun := max{l — £, gi, g2} < 1. This concludes the proof of fl4T]) . 

Step 4. Recall that ||h^+i||L°°(r2 ) < ||^o||l°°(o)- If C ||ho||L°°(o) < 1/4, estimate (ITT]) proves 

A^+i < giinA^- (1/4 - C ||/i£+i||^oc(f2)) |||fo+i -follP < TiinA^ for all i e Mq. 

With reliability fl22l) and osc^ < rjj from (ETj), induction on n proves 

7 h£+n < A^+n < gun Af < gHa(Crei + 7 + /^) vl all i,ne No. 

This proves linear convergence fl32l) with Cun = (Crei + 7 + /i)7~^- 
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Step 5. Together with the triangle inequality |||Mfc+i —< 2|||m —Mfc+i||p + 2|||-u — 
reliability fl2^ . and linear convergence ([32]), the geometric series yields 

^ ^+1 eij ^+1 Isa , 

'^\\\uk+i-Uk\f <4:^\\\u-Uk\f < ^ vI^qL^vI 

k=£ k=i k=l j=0 j=0 

This concludes the validity of the general quasi-orthogonality (A3). □ 


3.5. Auxiliary results. For the convenience of the reader, this section collects some 
well-known properties of the FVM which are exploited in the subsequent proofs. 

Lemma 13. With x* G 'P^O'x) being the eharacteristic funetion ofVi & T^, we define the 
interpolation operator 

i; : C(n) ^ V'‘(T;), := 5^ v{a,)x:. 

0,i EA/*x 

Then, for allT eTx, F E Ft, and E S^{Tx), it holds 

I {vx — F^Vx) dx = 0 = I {vx—F'^Vx) ds = 0, (42) 

J T J F 

Ikx ~ ||i, 2(7’) < hrll Vwx ||i,2(7’), (43) 

llwx -X>x||L2(ir) < Ch;^'^^||Vnx||L2(T)- (44) 

In particular, it holds F’^Vx E V^iTf) for all E Sq{Tx)- The constant C > 0 depends 
only on a-shape regularity of Tx ■ 

Proof. The proof of (l42l) is based on the construction of Ff and exploits that Vx is a piecewise 
linear function on 71<. A proof for (H3i) is found in |Erain[ Lemma 1.4.2], and (1441) follows 
from (1431) and the trace inequality. □ 


The following lemma is a key observation. For discrete ansatz and test spaces, it allows 
to understand the FVM bilinear form as a perturbation of the bilinear form of the weak 
formulation. The proof is given in |ELLn21 lChan2[ IEral2] for Lipschitz-continuous A, but 
transfers directly to the present situation, where A satishes (|7|)-lj8|). 

Lemma 14 f |ELL02l lCha021 lEral2] L It holds 

|Xl(nx, Wx) - Ax{vx,Fl,Wx)\ < Cbii ^ hr lll'^x |||t|||wx|||t for all Vx,Wx E Sq(Tx). (45) 

Ter 

Moreover, there exists some constant H > 0 such that Cbii||hx ||loo(o) < H implies 

Ax(Vx,Fy^Vx') P C*stab IWvxlW^ for allvx E Sq{Tx). (46) 

In particular, this proves that the FVM system (1T4]) has a unique solution Ux E Sq{Tx)- 
While H > Q depends only on the assumptions (I7j)-(IH|) on A, the constants Cbii and (Fstab 
depend additionally on a-shape regularity ofTx- □ 
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3.6. Proof of Theorem!^ The proof is split in several steps: 

Step 1. For arbitrary Vx,Wx G iSq( 7^) and ta* := X^Wx, we prove the identity 

A{vx,Wx)-Ax{vx,w*^) = f(divxAVt>x , w* -Wx)t - (AVwx ■ n,w*^-Wx)aT\r^ ■ (47) 

First, elementwise integration by parts for the bilinear form A.(nx,'?i'x) leads to 
A(ux,Wx)= ^ (AVwx, Vm;x)t = ( - (divxAVux , Wx)t + (AVwx ■ n,M;x)ar\r), 

TeTx TeTx 

since tCxIr = 0. Second, we rewrite the FVM bilinear form Ax(nx,'M^x)- Note that tc* does 
not jump across facets F E J^x- Therefore, 

Ax{vx,w*^) = '^2 ^xki / (—AVfx) • nds = — E E (AVux ■ n,wl,)TndVi- 

ai€Af^ TeTx aieAfrXr 


Note that A/t\F can be replaced by A/r, since tCxIr = 0. Integration by parts thus yields 
Ax(wx,w;) = -^ ^ (AVux ■ n,w;;)a(rny) - (AVux ■ n,w;;)ar\r) 

TeTx QiGA/y 

= ^ - (divxAVux , Wx)t + (AVux ■ n, w* jaTVr). 

TeTx 

The difference of the above estimates prove fl47p . 

Step 2. For arbitrary Vx,Wx E Sq{Tx) and tc* := X*tCx, h holds 

{f,Wx - Wx) - {Ax{vx,Wx) - A{vx,Wx)) < C'oscx(fx) ||Vwx ||l2(o), (48) 

where C > 0 depends only on a-shape regularity of 77<: With fl47)) . the definition of the facet 
residual flT5|) . the L^-orthogonalities (H2|) and the Cauchy-Schwarz inequality, we see 

(/,«;* - Wx)n - {Ax{vx,wl) - A{vx,Wx)) 

= ^ ((1 - nx)(/+ diVxAVfx) , W* - Wx)t - ((1 - nx)Jx(t’x), W* - Wx)f 

TeTx FeT^ 


< OSCx(Vx) ( 


Wx 


Wy 


ih(T)) 


1/2 


+ oscx(nx) ( 


hj}\ 


Wx 


Wy 


|2 

\L^{dT\V) J 


TeTx TeTx 

With fl45]) - (HT)) . we conclude fl45]) . 

Step 3. The FVM solution Ux satisfies 

|||mx — Wx III < |||m — fx III + oscx(vx) for all fx e >So(7)<), (49) 

where C > 0 depends only on a-shape regularity of Tx and the assumptions (I7])-([H]) on A: 
With u being the weak solution, we first note the identities 

0 = {f,Wx)n - A{u,Wx) = [if,Wx)n - A(r’x,tWx)] - A{u - Vx,Wx) for all Vx E Sq{Tx). 
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For sufficiently fine Tx, Lemma [14] applies. Choose Wx ■= Ux — Vx and ta* := I^Wx- Then, 


diet 


Ux - Vy 


r ^ -^xiux - Vx,wl,) 


Gsli 


= [(/,'iax)n - ^x(wx, w*)] - [{f,Wx)n-A{vx,Wx)] + A{u - Vx,Wx)- 

Combining this with fl48|) and norm equivalence ||VtCx||L2(r2) — |||'iax|||, we obtain 


\Ux - Wxlf < OSCx(nx) llkxlll + \\\u - Vy. 


\w 


xMh 


where the hidden constant depends only on a-shape regularity of Tx and the assumptions (j7])- 
([H]) on A. By choice of tCx, we conclude flT^ . 

Step 4. Let Ux £ SliTx)- We employ (Bl’) with % = Tx and = u^- Combining this 
with the triangle inequality and (1491) . we see 

l l49ll 

|||m-Mx||| +0SCx < lllu-nxlll + |||mx -^^x||| +osCx(nx) < |||M-nx||| +osCx(nx)- 

Altogether, this proves 

|||m — Mx III + oscx ^ min (|||n — Uxlll + oscxfux)) < |||w — Mx||| + oscx- 

Vx&S^iTx) 

Reliability (1221) and efficiency (125]) together with (12T]) imply px — lll'*^~'Wx|||+ oscx- This 
concludes (124|) . For the equivalence 


\u — u 


FEMi 

X I 


+ OSCx (U 


FEM^ 


X V“x 


~ mm 


Vx (Tx) 

the reader is referred to |FFP141 Lemma 5.1]. This also concludes 


U-Vx\\\ +OSCx{Vx)), 


□ 


3.7. Proof of Theorem [3l Let Vx G iSg(7^). For T E Tx, the triangle inequality shows 

osCx(T,nx)^ < 11(1 - Llx)/||i2(7^) + ||(1 - nx)divxAVux|li2(7^) 

+ hx 11(1 — Llx)</x(nx)||i2(a7’\r). 

With Step 1-2 from the proof of Lemma [TUI we thus see 

OSCx (Ux) ^ II ^x (1 ~ nx)/|lL2(0) + ||hx Vux |Il2(q). 

With ([24]) . we obtain 

ll24ll 


m — Mx + OSCx — 


mm 

V X £* 5(1 (Tx ) 


\u - Uxlll + oscx(nx)) 


< 


||hx(l - nx)/||L 2 (Q) + min (|||m - UxIII + ||hxVux|| l2(o))• 

Vx€S^{Tx) 

This proves ([271) . In particular, norm equivalence |||m — UxlH — ||V(m — nx)||T2(f7) implies 

Illn-Wxlll + 0 SCx < ||hx|U°o(o)(||/||T2(Q) + ||VM||L2(f^)) + min |||M-nx|||- 

VxGS^iTx) 

From this, we also conclude ([28l) - ([29l) . 


□ 


3.8. Discrete reliability. The main result of this section is the following variant of the 
discrete reliability (A4). 
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Proposition 15. Let Tx G refine(7^) he an arbitrary refinement of To & refine(7o) and 
suppose that the corresponding discrete solutions Ux or Uo exist. Then, 

|||mx - MoIlP < Cu\ ^ h%\\\u-x - Uo\\\t + Cdir (50) 

TgTx T&TZo 

where IZo := {T G 7^ : 3T' G To\Tx T f] T' 7 ^ 0}, consists of all refined elements 
To\Tx plus one additional layer of neighboring elements. In particular, the discrete reliabil¬ 
ity (A4) follows provided thatTx is sufficiently fine, i.e., Cbii||/ix |lioo(Q) < 1 / 2 . The constants 
Cbii, Cdir > 0 depend only on 12, the assumptions 0 -l|8|) on A, and on a-shape regularity of 

To. 

The proof of Proposition (TS] relies on two properties of the volume and facet residual, i.e., 
an orthogonality property flSTl) and a discrete defect identity fl5^ of the FVM bilinear form. 

Lemma 16. LetTx G refine(7^) be an arbitrary refinement of To G refine(7o) and suppose 
that the corresponding discrete solutions Ux or Uo exist. Then, there holds 

J2{Jo{uo),v:)f = 0 forallvleV^.iT:) (51) 

as well as 


^ (i?o(Mo) ,'yx)T - ^ (7 o(mo) ,r’x)F = Ax(mx - Mo, Wx) for all v*^ & . (52) 




FeJf 


Proof. The proof of fl5T|) is well-known and found, e.g., in |CLT051 [Eral0| lEral3] . The proof 
of fl5^ is adopted from [ZoulO] for an arbitrary refinement Tx G ref ine(7^); The divergence 
theorem shows for all boxes V G 75* that 


y / divxAVMo(ir= y. / Jx{uo)ds-\- / AVwo-nds. 

T'gTx ^ JdV 


(53) 


Let u* G VoiTf). We multiply the above equation by and sum over all V ^ Tf. With 
divxAVwo = div^AVtio, the left-hand side then reads 


V'GT; T'gTx 


iT'nV 


divxAVtiodT = (divxAVMo,M*)Q = ^ (div^AVwo,T x)t- (54) 


TgTo 


Since AVmo is continuous in T G 75 and Jx(mo) = Jo{uo) on F G T^, it holds 


f Jx{Uo)dS= ^ (J; 

V'&T* C'eT^,_„ F'&T^ 


x[Uo),Vx)f' — ^ ^ iJo[Uo) ,Vy^)F. 
T'gT“ TgT^ 


By dehnition (IT^ of Ax(-, •)> l^e identity (|53l) becomes with (j54|) and (|55l) 

(divoAV mo ,v*x)t= (Jo(Mo) ,Mx)f - Ax(Mo,Mx). 

T£To FgTP 

On the other hand the FVM formulation (IT^ yields 

(/,Mx)o = Ax(Mx,Mx). 

Adding fl56l) - fl5711 . we conclude the proof. 
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(55) 


(56) 


(57) 

□ 








The following Poincare- and trace-type inequalities play a key role to estimate quantities 
over the elements of the dual grid. 

Lemma 17. For each box Vi E let Oi G Afx be the corresponding node. Define 

HI ■. L\a) ^ vHt;), (n*t,)|v, 

[0, ifOiEM^. 

Let V eT^ and ( E Dv,x- Then, there holds, for all v E Hq{Q), 

llv - n>||i 2 (y) < Cdiam(l/) ||Vn||L 2 (y), (58) 

||n - n* n|| 2 , 2 (^) < Cdiam(l/)^/^ ||Vn||L 2 (y). (59) 

The constant C > 0 depends only on the a-shape regularity of Tx ■ 

Proof. The set Tx\v ■= {V FT : T E Tx with V P\ T 7 ^ 0} is a partition of V into 
quadrilaterals in 2D and cuboids in 3D, respectively. In 2D each quadrilateral can itself be 
divided into two triangles. In 3D each cuboid can be divided into three pyramids (with the 
center of gravity of T as top). Note that a quadrilateral ( E J-'v,x builds the base of one 
pyramid. This gives rise to a triangulation Zy^x of V] see Figure [Hand Figure [2] for 2D and 
3D, respectively. 

Choose Z E Zy^x with C, C dZ. Note that Zy^x is a'-shape regular, where a' depends only 
on cr, and that the box V is just the node patch of the corresponding node a E Mx with 
respect to Zyx. If a G Nff, let vz '.= {1/\Z\) J^vdx denote the piecewise integral mean. If 
a e Ml we define := 0, since then D fl F has positive measure. In either case, it holds 

lln - n>||L 2 (y) < ||n - Vz\\L^(y) < diam(Z) || Vn||L 2 (y), 

where the hidden constant depends only on a' and hence on a; see |DS80] . With diam(Z) < 
diam(I/), the Poincare-type inequality (158|) follows. 

The trace inequality, a scaling argument, and diam(I/) ~ diam(Z) lead to 

lblU 2 (c) < diam(Z)"^/ 2 ||n||i 2 ( 2 ) -h diam(Z)^/ 2 || Vn||L 2 (^) 

< diam(I/)"^/^|ln||i2(y) -F diam(I/)^/^|| Vn||L2(y). 

Combining this with the Poincare-type inequality fl5S]) . we obtain 

ll» - n,*K||i2K, < diam(V')‘''"||Vt,|U.,v,. 


This concludes the proof. 


□ 


Proof of Proposition IT5L To abbreviate notation, let := Ro{u,f) and Jo := Jo(mo) 
denote volume residual (ITTD and facet residual (ITS]) with respect to the discrete solution u^. 
For arbitrary Vx E Sq(Tx), ^ '^oO'x), and v* E Vq{%), (IS^ and (1ST]) of Lemma ITS] 
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(a) Mesh To and Tx (b) Dual mesh 7^*. 
(inch dashed lines). 



(c) Dual mesh 7^* . 


(d) Set TZo- 


Figure 4. In (a) we see the coarse mesh To for 2D. The dashed lines show 
the rehnement and build the rehned mesh Tx- In (b) and (c) (gray boxes) we 
see the dual mesh of the rehned areas notated by 7?.* and 7?.*, respectively. 
Finally (d)| shows the elements T E % which build the set TZo ■= {T E Ua\a C 
diTo\Tx)} in this example. 


together with the mesh relation fllip show 
Ro) fx) 

= A{Ux - Uo,Vx) - Ax{Ux - Uo,V*^) + ^ (i?o,nx - K)t - X] (4,^* - vI)f 

t&T f&jT 

= A{u^-Uo,v,,)-Ax{ux-Uo,vD+ ^ {{Ro,v*^-vl)v - ^ (4,nx-Oc\r); 

(60) 

see Section [2751 Figure and Figure for the dehnition of Tv,o- 

Next, we note that the discrete ansatz spaces are nested, while the discrete test spaces 
are not. However, in the non-rehned area TnTx the shape of the dual grid elements is the 
same. We use this to truncate the sum of fl60ll . To get the hnal sum over TZo in fl50|) . we 
have to dehne the functions w * and v* appropriately to apply Lemma [13] and Lemma [T71 
respectively. To formalize this, we dehne TZl := To\T^ and TZ*^ := T^\TT i-^-, the dual 
mesh of the rehned areas; see Figure 01 for a 2D illustration. Note that 

U U u (61) 

v&n% u'eR* 

Consider the transition area TZo\{To\Tx) = {T E T^Tx '■ 3T' E To\Tx T fl T' 7 ^ 0} 
(second row of triangles in Figure |3|) which consists of all non-rehned neighbors of a rehned 
element. For all T E TZo\(To\Tx), it holds 

{FnT : F e 7^;} = {FnT : V eTZ*^}, 
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i.e., the shape oiVE Til coincides with the shape of some V G Tl*^ in the transition area. 
Let fx := Uy, — ^ Sl{Ty). Choose n* := X^Vy G V^iTy). Dehne f* G V^iT^) by 


(n:nx)iv ifCG7^:, 
otherwise. 


i.e., v: 


= vl within the white area 


For V G T^yJl l = %* r\ Ty, this implies vUv = vl\v, 
of Figure l ^b) and 0|( 
replace 7^* by 7?.*. Together with fl45|) from Lemma [14] for the bilinear forms, we get 


We use this observation to truncate the sum over 7^* 


m 


and 


A{u^ 


Uo.V 


O’) J ^ 


< Cbii ^ hj 




Uc 


T T> 


TgTx 


+ 5^ [{R,,vi-v:)v- J2 

V C^^V,o 



(62) 


Next, we estimate the sum over T G 7^ by the Cauchy-Schwarz inequality. Furthermore, we 
add Vy — Vy and use fIM]) to rewrite the sum over the boxes V G Til (ES]): 


.4(mx -u^,vy) < ( X] ^tIII^x -MoIIIt) III^^xIII 

^TgTx 

+ ((^^Wx-Ov- (-^^Wx-Oar) 

V £71% 

+ X] [{R<,,vl-Vy)v' - ^ (4Wx-'f^x)c'\r)- 

V'£ni C'eWCx 

C'cFej=-o 


(63) 


Note that 7^v',x contains also parts of facets from Xy which are not needed here and which 
are avoided by C F G F*. To abbreviate notation, let hy '■= diam(l/) and note that 
(T-shape regularity implies hy — hr for all V E and T E % with V O T 7 ^ 0. Next, 

we estimate the two sums over Til F*: First, with (158|) and fl5^ of Lemma [T7I and 
K\v = n^'^x |v for all V E Til, ^be Cauchy-Schwarz inequality yields 

(^{R^,Vy -vl)y - Y (4Wx-Oc\r) 


< 

rv_/ 


(E hi\mhmy'R{ E E 

- V£'R.% V£'R.%C£j^v,o 


LHCX) 


1 / 2 ' 


E ii^*' 


vent 


||2 \ 
x|lL2(y) j 


1/2 


With UvGTes ^ Utgt^o hence obtain 


< 

rs_/ 


Y^ (br ll-RollL2(r) + b-T ||bo||i2(9-r\r)) 


11/2 


'-TGTeo 


TeTio 


(64) 
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Note that U^/'G 7 ^* ^ UreTeo Then, with (143|) and (jHj) of Lemma fT3l and n* = XyV 

we get as before 




.^X - ^xjy' 


(Jo,<-'yx)c 


'\r 


v’ani 


C'cFeJ='o 


< 

rs_/ 


i^T ll-Ro|lL2('r) + hr ||To|1^2(g2’\r)) 


'-Te'Ro 


1/2 


Vx = 




1/2 


tgt^o 


Combining with flB^ . we obtain 


j\.(Ux U^, Vx) 


< 


S'* 


2 Ml Ill2 ^ 

rp I I I W I I I X' 


1/2 


J 




TgTx ^G7^o 

Finally, ellipticity of A{-, ■) and the choice of Vx = Uy — show 


1 / 2 ' 


Ux - Ua 


< 


A{Ux Uo) Vx) 


< 


^ ^ hj’ |||mx HoIIIt “I" ^ ^ ho(T') 


'-TgTx 


Te7^o 


1/2 


X "^X ) 


(65) 


-Ux - Mo 


This proves fl50|) and concludes the proof. 


□ 


3.9. Proof of Theorem O Suppose that the initial triangulation To is sufhciently hne 
such that the following assumptions (i)-(iii) are satished: 

(i) For all Tx G ref ine(7o), the FVM system (IT^ is well-posed. In particular. Lemma [9] 
proves that stability (Al) and reduction (A2) are satished. 

(ii) Proposition [TT] is valid and, in particular, the general quasi-orthogonality (A3) is 
satished. 

(iii) The constant Cbu from Proposition ITS] satishes Cbii||ho||Loo(o) < 1/2, so that Propo¬ 
sition [TSl in fact, proves the discrete reliability (A4). 

Finally, let A4.i C 7^ be a set of minimal cardinality which satishes the Dorher mark¬ 
ing (l30|) for the error estimator. Then, the additional assumption of Theorem [7] (ii) and 
the choice of the marked elements A4/ C C 7^ in Algorithm 0] imply that < 

Cmns < C'MNsC'mark #A4£. Altogether, the assumptions of |CFPP14l Theorem 4.1] are 
fulhlled, and (157]) -(I55 ]) follow for our adaptive FVM of Algorithm 01 □ 

4. Numerical experiments 

In this section, we illustrate the performance of Algorithm 0] with 9 = 0.5 = 9' for two 
examples. In extension of our theory, we consider the model problem ([6]) with inhomogeneous 
Dirichlet boundary conditions. The numerical experiments are conducted in Matlab on a 
standard laptop with a dual core 2.8 GHz processor and 16 GB memory. 

4.1. Experiment with smooth solution. On the square O = (—1,1)^, we prescribe 
the exact solution u{xi,X 2 ) = (1 — lOx^ — with x = {xi,X 2 ) G We choose 

the diffusion matrix 

10-l-cosxi 9 x 1 X 2 

9xiX2 10 -|-sinx 2 
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(a) To (16 elements). 



(c) Tie (3838 elements). 



(b) Ti (216 elements). 



x-i 

(d) Solution (Tie). 



Figure 5. Experiment with smooth solution from Section IFTI Initial trian¬ 
gulation To with NVB reference edges as well as adaptively generated meshes 
Ts resp. Tie, and discrete FVM solution calculated on Tie- 


so that ([7]) holds with Amin = 0.82293 and Amax = 10.84096. The right-hand side / is calcu¬ 
lated appropriately. The uniform initial mesh consists of 16 triangles; see Figure!^ 


In Figured 1(b) and d 1(c) we see an adaptively generated mesh after 8 and 16 rehnements, 
respectively. Figured 1(d) plots the smooth solution on the mesh Tie- As u is smooth, uniform 
and adaptive mesh-rehnement lead to the optimal convergence order CT(iV“^A) with respect 
to the number N of elements; see Figure [61 The oscillations are of higher order and decrease 
Table 


with 0 {N-^). 
Theorem [7] (ii 


shows the experimental validation of the additional assumption in 


that marking for the data oscillations is negligible. 


4.2. Experiment with generic singularity. On the L-shaped domain O = 
(-l,in(|0.1]x[-1.0]) , we prescribe the exact solution u{xi,X 2 ) = sm.{2Lp/?t) in polar 

coordinates r G Mq, (p G [0,27r[, and {xi,X 2 ) = r(cos sin <yc). Then, u has a generic singu¬ 
larity at the reentrant corner (0, 0), which leads to u G iT^+^A-£()y2) for all e > 0. We choose 
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Figure 6. Experiment with smooth solution from Section IFTI Error in the 
energy norm Ei := |||m — m^IH, weighted-residual error estimator rji, and data 
oscillations osc^ for uniform and adaptive mesh-rehnement. 


the diffusion matrix 


A = 


5 + cosxi 

{xl + xD^ b + {x\ +x\)si\iX 2 


xlf 


so that ([7]) holds with Amin = 0.46689 and Amax = 5.14751. The right-hand side / is calculated 
appropriately. The uniform initial mesh consists of 12 triangles. Some further adaptively 
generated meshes together with a plot of the discrete solution are shown in Figure [71 

For uniform mesh rehnement, we observe the expected suboptimal convergence order 
of while adaptive mesh-refinement regains the optimal convergence order of 

see Figure |H1 As in the experiment of Section 14.11 the oscillations are of higher 
order 0{N~"^). See Table Ql^b) for the experimental validation of the additional assumption 
in Theorem [7| (ii) that marking for the data oscillations is negligible. 
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energy norm Ei := |||m — m^IH, weighted-residual error estimator rji, and data 
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(a) Section ITT] 



(b) Section IT? 


i 

#Ti #Me/#M'; 

OSc{M'l,Y / oscj 

£ 

WtTTmJWMT 

osc{MJY/oscj 

0 

16 

1.000 

0.634 

0 

12 

1.667 

0.143 

1 

22 

1.000 

0.613 

1 

18 

1.750 

0.115 

2 

28 

1.000 

0.704 

2 

26 

1.400 

0.108 

3 

32 

1.000 

0.769 

3 

35 

1.222 

0.062 

4 

40 

1.214 

0.338 

4 

56 

1.200 

0.104 

5 

78 

1.111 

0.449 

5 

78 

1.643 

0.028 

6 

112 

1.133 

0.292 

6 

no 

1.350 

0.135 

7 

156 

1.119 

0.410 

7 

148 

1.161 

0.290 

8 

216 

1.062 

0.393 

8 

204 

1.111 

0.268 

9 

331 

1.198 

0.263 

9 

274 

1.048 

0.423 

10 

460 

1.014 

0.474 

10 

370 

1.168 

0.223 

11 

660 

1.049 

0.371 

11 

525 

1.069 

0.324 

12 

944 

1.027 

0.430 

12 

704 

1.063 

0.296 

13 

1,340 

1.025 

0.404 

13 

961 

1.015 

0.442 

14 

1,914 

1.019 

0.383 

14 

1,314 

1.003 

0.475 

15 

2,752 

1.026 

0.374 

15 

1,784 

1.037 

0.345 

16 

3,838 

1.015 

0.358 

16 

2,451 

1.000 

0.639 

17 

5,428 

1.003 

0.449 

17 

3,305 

1.015 

0.417 

18 

7,430 

1.013 

0.359 

18 

4,562 

1.000 

0.595 

19 

10,572 

1.003 

0.445 

19 

6,161 

1.001 

0.482 

20 

14,462 

1.019 

0.322 

20 

8,344 

1.011 

0.440 

21 

20,264 

1.004 

0.431 

21 

11,316 

1.000 

0.635 

22 

27,532 

1.004 

0.455 

22 

15,249 

1.000 

0.528 

23 

38,402 

1.010 

0.323 

23 

20,631 

1.000 

0.577 

24 

52,366 

1.000 

0.539 

24 

27,742 

1.014 

0.451 

25 

72,386 

1.007 

0.401 

25 

37,566 

1.000 

0.655 

26 

98,144 

1.000 

0.509 

26 

50,139 

1.011 

0.437 

27 

135,076 

1.004 

0.445 

27 

67,722 

1.000 

0.571 

28 

184,006 

1.000 

0.605 

28 

90,543 

1.000 

0.523 

29 

251,668 

1.002 

0.475 

29 

121,136 

1.005 

0.471 

30 

341,940 

1.001 

0.488 

30 

163,221 

1.000 

0.715 

31 

461,354 

1.000 

0.616 

31 

216,681 

1.025 

0.361 

32 

634,922 

1.004 

0.415 

32 

292,527 

1.000 

0.545 

33 

852,264 

1.000 

0.663 

33 

389,411 

1.000 

0.582 

34 

1,171,426 

1.002 

0.465 

34 

521,975 

1.013 

0.437 

35 

1,567,542 

1.000 

0.611 

35 

699,195 

1.000 

0.678 

36 

2,150,232 

1.000 

0.521 

36 

928,417 

1.012 

0.418 

37 

2,893,626 

1.000 

0.652 

37 

1,246,972 

1.000 

0.561 

38 

3,932,562 

1.000 

0.593 

38 

1,658,877 

1.000 

0.585 

39 

5,335,740 

1.000 

0.493 

39 

2,224,754 

1.003 

0.481 

40 

2,959,035 

1.000 

0.659 


Table 1. Experimental results on marking strategy: We compute Cmns •= 
< 2 and see that the additional assumption in Theorem [7] (ii) is 
experimentally verihed. In addition, we compute 9' := osc£(Wl^)^/osc^ > 0.02, 
i.e., the choice 9 = 0.5, 9' = 0.02 would guarantee Aii = Ai] in Algorithm 01 
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